Latest inflation model constraints from cosmic microwave background measurements 
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We present an update of cosmological constraints on single-field inflation in light of the Wilkinson 
Microwave Anisotropy Probe satellite mission five-year results ( WMAP5) . We find that the cosmic 
microwave background data are quite consistent with a Harrison- Zel'dovich primordial spectrum 
\^ . with no running and zero tensor amplitude. We find that the three main conclusions of our analysis 

^ ■ of the WMAP three-year data (WMAP3) are consistent with the WMAP5 data: (1) the Harrison- 

Zel'dovich model is within the 95% confidence level contours; (2) there is no evidence for running of 
I I , the spectral index of scalar perturbations; (3) From the WMAP 5 data alone, potentials of the form 

i-^ ' V oc (f)^ are consistent with the data for p — 2, and are ruled out for p = 4. Furthermore, consistent 

with our WMAP3 analysis, we find no evidence for primordial tensor perturbations, this time with 
Q I a 95% confidence upper limit of r < 0.4 for the WMAP5 data alone, and r < 0.35 for the WMAP5 

■ data taken in combination with the Arcminute Cosmology Bolometer Array (ACBAR). 

, PACS numbers: 98.80.Cq 
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' Inflation ,1] is the dominant paradigm for understanding the initial conditions for structure formation and for 
Cosmic Microwave Background (CMB) temperature anisotropics. In the inflationary scenario, primordial density 
and gravitational-wave fluctuations are created from quantum fluctuations, "redshifted" beyond the horizon during 
I an early period of superluminal expansion of the universe, then "frozen" [3, S fl- Perturbations at the surface of 
. last scattering are observable as temperature anisotropies in the CMB, as first detected by the Cosmic Background 
L| ' Explorer satellite. The latest confirmation of the inflationary paradigm has been recently provided by the five-year 
• !-( data from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite [E H, 0, B S [111 ■ The WMAP collaboration 
has produced new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz based on the 
JH . first five years of the WMAP sky survey. The five-year maps incorporate several improvements in data processing made 
■ - - ' possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight 
beam response. WMAP data support the inflationary model as the mechanism for the generation of super-horizon 
curvature fluctuations. 

The goal of this paper is to make use of the recent WMAP flve-year data to discriminate among the various single- 
field infiationary models. As such, this paper represents an update of our pre vious analyses [Tll.fl^ of the first-year 
and three- year WMAP data, respectively. The reader is referred to Ref. for a more detailed discussion of the 
inflationary model space, and for more extensive references. 

We compare our results to those reported by the WMAP team in Komatsu et al., Q and in Dunkley et al. Q- As 
first pointed out in Ref. for single-field infiation models, the relevant parameter space for distinguishing among 
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models is defined by the scalar spectral index n, the ratio of tensor to scalar fluctuations r, and the running of the 
scalar spectral index dn/d\nk. The Harrison-Zeldovich (H-Z) point {n = l,r = 0) is comfortably within the 95% 
C.L. contour of Ref. [1], and just at the borderline of the 95% C.L. contour of Ref. [3]. This small difference appears 
to be statistically insignificant. In our analysis using the WMAP5 data with slightly different priors than used by the 
WMAP team (discussed in Sec. Ill) we find that marginalizing over all parameters except n and r in the WMAP5 
data, the likelihood of the H-Z point r = 0, n = 1 is 0.0988, corresponding to exclusion of the H-Z model with 
confidence level of 93.8%. Marginalizing over all parameters except n {i.e., the single-parameter confidence limits 
on n for the same parameter set), the likelihood of the point n = 1 is 0.4226, corresponding to exclusion of the 
scale-invariant spectrum with a confidence level of 81.1%. These figures are highly dependent on our choice of priors: 
if we use a six-parameter fit with a prior of r = 0, the likelihood of the point ti = 1 is 0.07392, corresponding to the 
exclusion of the scale-invariant spectrum with a confidence level of 97.8%. (In all cases, the maximum likelihood is 
normalized to unity.) Our conclusion is that the H-Z model is not preferred, but neither is it strongly disfavored. 

The paper is organized as follows: In Sec. II we will quickly review single-field infiation models and their observables. 
In Sec. Ill we define the inflationary model space as a function of the slow-roll parameters e and rj. In Sec. IV we 
describe our analysis method as well as our results. Since a study of the implications of the WMAP5 data for single 
field models of inflation has been already performed by the WMAP collaboration themselves |5(| , we will also specify 
the differences between our analysis and theirs. In Sec. V we present our conclusions. 



II. SINGLE-FIELD INFLATION AND THE INFLATIONARY OBSERVABLES 



Inflation not only explains the large-scale homogeneity of the universe, but also provide a mechanism for explaining 
the observed level of inhomogeneity as well. During inflation, quantum fluctuations on small scales are quickly 
redshifted to scales much larger than the horizon size, where they are "frozen" as perturbations in the background 
metric. The metric perturbations created during inflation are of two types: scalar, or curvature perturbations, which 
couple to the stress-energy of matter in the universe and form the "seeds" for structure formation, and tensor, or 
gravitational-wave perturbations, which do not couple to matter. Both scalar and tensor perturbations contribute 
to CMB anisotropics. Scalar fluctuations can also be interpreted as fluctuations in the density of the matter in the 
universe. Scalar fluctuations can be quantitatively characterized by the comoving curvature perturbation P-ji. As 
long as slow-roll is attained, the curvature (scalar) perturbation at horizon crossing can be shown to be Jj] 
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Here, denotes the slow-rolling scalar field dominating the energy density of the Universe during infiation, the so- 
called inflaton, H is the Hubble rate and mpi = 1.2 x 10^^ GeV is the Planck scale. The slow roll approximation is 
consistent if both the slope and curvature of the inflaton potential V{(j)) are small (in Planckian units), V , V" ^ V. 
In this case the slow-roll parameter e can be expressed in terms of the potential as 
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We will also need in the following a second "slow-roll parameter" 77 defined by 
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Slow roll is then a consistent approximation for e, 77 ^ 1. 

The fluctuation power spectrum is, in general, a function of wavenumber fc, and is evaluated when a given mode 
crosses outside the horizon during inflation, k = aH. Outside the horizon, modes do not evolve, so the amplitude 
of the mode when it crosses back inside the horizon during a later radiation- or matter-dominated epoch is just its 
value when it left the horizon during inflation. Instead of specifying the fluctuation amplitude directly as a function 
of k, it is convenient to specify it as a function of the number of e-folds N before the end of inflation at which a mode 
crossed outside the horizon. 

The scalar spectral index n for P-ji is defined by 
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so that a scale-invariant spectrum, in which modes have constant amplitude at horizon crossing, is characterized by 
n = 1. 

The power spectrum of tensor fluctuation modes and the corresponding tensor spectral index is given by [l| 
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The ratio of tensor-to-scalar modes is then Pt/P-r — 16e, so that tensor modes are negligible for e <C 1. In the 
limit of slow roll, the spectral indices n and ut vary slowly or not at all with scale. We can write the spectral indices 
n and tlt to lowest order in terms of the slow roll parameters e and t] as 
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The tensor/scalar ratio is frequently expressed as a quantity r, which is conventionally normalized as 

Pt 
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The tensor spectral index is not an independent parameter, but is proportional to the tensor/scalar ratio, given to 
lowest order in slow roll by ut — — 2e = — r/8. This is known as the consistency relation for inflation. A given 
inflation model can therefore be described to lowest order in slow roll by three independent parameters: P-jz, Pt, and 
n. If we wish to include higher-order effects, we would have to include a fourth parameter describing the running of 
the scalar spectral index, dn/dlnk. 



III. ANALYSIS AND RESULTS 



The method we adopt is based on the publicly available Markov Chain Monte Carlo (MCMC) package cosmomc 
[3]. We sample the following eight-dimensional set of cosmological parameters, adopting flat priors on them: the 
physical baryon and CDM densities, ujb — fl^h^ and ujc = f^c^^, the ratio of the sound horizon to the angular diameter 
distance at decoupling, 9s, the scalar spectral index, its running, and the overall normalization of the spectrum, n, 
dn/dlnk and A at some pivot scale k (either k = 0.002 Mpc~^ or fc = 0.017 Mpc^^), the tensor contribution r, and, 
finally, the optical depth to reionization, r. Furthermore, we consider purely adiabatic initial conditions, we impose 
spatial fiatness, and we use the inflation consistency relation to fix the value of the tensor spectral index ut- We also 
restrict our analysis to the case of three massless neutrino families; introducing a neutrino mass in agreement with 
current neutrino oscillation data does not change our results in a significant way. 

We include the five-year data |5|] (temperature and polarization) with the routine for computing the likelihood 
supplied by the WMAP team and available at the LAMBDA web site prttp:77iambda.gsfc.nasa. gov/) We also perform an 
analysis based on the WMAP five-year data in combination with results from the Arcminute Cosmology Bolometer 
Array ( ACBAR) [l^ . The MCMC convergence diagnostic is done on 8 chains though the Gelman and Rubin "variance 
of chain mean" /"mean of chain variances" R statistic for each parameter. Our ID and 2D constraints are obtained 
after marginalization over the remaining "nuisance" parameters, again using the programs included in the cosmomc 
package. Furthermore, we make use of the HST measurement of the Hubble parameter Hq — 100ft- km s~^Mpc^^ 
[iSj by multiplying the likelihood by a Gaussian likelihood function centered around h = 0.72 and with a standard 
deviation a — 0.08. Finally, we include a top-hat prior on the age of the universe: 10 < to < 20 Gyrs. We perform 
the analysis in two ways: first, we constrain parameters with a prior assumption of dn/dlnk — 0, i.e., an exact 
power-law perturbation spectrum, and second, with dn/dlnk a free parameter in the MCMC. We include lensing in 
the calculation, and marginalize over the Sunyaev-Zel'dovich amplitude using the WMAP V-band template. 

All fits are performed at a pivot scale of fc = 0.017, which is different from our previous analysis and the analysis 
of Dunkley et al. and Komatsu et ai, all of which are normalized to a pivot scale of A: = 0.002. The reason for the 
different choice of pivot scale is to mitigate the parameter degeneracy between n and dn/dlnk [l7|, which significantly 
improves convergence of the MCMC code. In the case of dn/dlnk — 0, the choice of pivot scale has no effect on the 
spectral index n, since a change of pivot affects only the normalization: 
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However, the change of pivot results in a rescaling of the tensor/scalar ratio r which is dependent on r and n, since 
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and therefore, using the consistency relation ~ 

Figure [1] shows 68% and 95% confidence limits on r and n in the case of a no-running prior, dn/d\nk = 0. Clearly 
the choice of the pivot scale makes very little difference in this case (and since it is just a normalization for r, no 
difference at all in the r = case.) 

In Fig. [2] we show the 68% and 95% confidence limits on r and n in the case where we allow running of the spectrum 
dn/d In fc as a free parameter. Again the WMAP5 alone data correspond to the open contours and WMAP5 + ACBAR 
to the filled contours. 

Finally in Fig. [3] we show the 68% and 95% confidence limits on dn/dhik as a function of 7i and r as a function of 
dn/d\nk. 



IV. CONCLUSIONS 



In this paper we presented an analysis of the recent WMAP five-year data set with an emphasis on parameters 
relevant for distinguishing among the various possible models for inflation. This analysis is an update to our similar 
analysis of the WMAP three- year data set In that analysis, we reached three main conclusions, which are either 
unchanged or slightly strengthened by the new data: (1) the Harrison-Zel'dovich model is within the 95% confidence 
level contour. We find that marginalizing over all parameters except n and r in the WMAP5 data, the likelihood of 
the H-Z point r = 0, n = 1 is 0.0988, corresponding to exclusion of the H-Z model with confidence level of 93.8%. 
Marginalizing over all parameters except n {i.e., the single-parameter confidence limits on n for the same parameter 
set), the likelihood of the point n = 1 is 0.4226, corresponding to exclusion of the scale-invariant spectrum with a 
confidence level of 81.1%. These figures are highly dependent on our choice of priors: if we use a six-parameter fit 
with a prior of r = 0, the likelihood of the point 71 = 1 is 0.07392, corresponding to the exclusion of the scale-invariant 
spectrum with a confidence level of 97.8%. (In all cases, the likelihood is calculated as the marginalized posterior 
probability, with the best-fit point normalized to unity. The confidence level is calculated from the likelihood assuming 
a Gaussian error distribution.) (2) There is no evidence for running of the spectral index of scalar perturbations. 
Figure [3] shows likelihood contours for dn/dlnk relative to the parameters n and r, calculated at a pivot scale of 
k = 0.017 to minimize parameter degeneracies, showing that the case of a pure power-law spectrum dn/dlnk = 
is fully consistent with all data sets. (3) Potentials of the form V (x (/>'' are consistent with the data for p — 2, and 
are ruled out by the WMAP five-year data alone for p = 4. This is a stronger conclusion that was possible with the 
WMAP three-year data, which was marginally consistent with at a 95% C.L. We find no evidence for a nonzero 
tensor/scalar ratio, with a 95% C.L. upper limit of r < 0.4 for the WMAP five-year data alone, and r < 0.35 for 
WMAP in combination with ACBAR. Our results are consistent with the results of Peiris and Easther in Ref. [l3| ■ 

Finally, we illustrate the consistency of CMB determinations of n and r over the years in Fig. [H The pre- WMAP 
status is well described by the combination of the COBE [19], BOOMERANG MAXIMA 21], DASI [22], and 
CBI [23] data sets. This pre- WMAP analysis includes also a prior on the optical depth r < 0.2. As illustrated by the 
figure, higher precision data prefers n < 1, but n = 1 is not strongly disfavored. 
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FIG. 1: 68% C.L. and 95% C.L. contours in the n, r parameter space for WMAP5 alone (open contours) and WMAP5 + 
ACBAR (filled contours), with a prior of dn/dln fc = 0. In the left figure the line segments show the predictions for V{cj)) = m'^cj)'^ 
(lower line segment) and V{(f>) ~ \(f>'^ (upper line segment) for the number A*' of e-folds before the end of inflation at which a 
mode crossed outside the horizon in the range A'' = [46, 60]. Also in the left figure the dashed lines show the 68% C.L. and 95% 
C.L. contours from our previous analysis of the WMAP 3-year data. The figure on the left is for a pivot scale of fc = 0.017 and 
the figure on the right is the same contours translated to a pivot scale of fc = 0.002 by Eq. (|10|l . The dotted line in the right 
panel is the result of a MCMC analysis run directly at a pivot of = 0.002 using the WMAP5 data alone. 
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FIG. 2: 68% and 95% confidence limits on r and n for WMAP5 alone (open contours) and WMAP5 + ACBAR (filled contours) 
allowing the possibility of a running spectral index. 
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FIG. 3; 68% and 95% confidence limits for dn/d\nk vs. n and r vs. dn/d\nk. The different contours are for different data sets 
as in Figs. 1 and 2. 
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FIG. 4: 68% and 95% confidence limits on r and n with a prior of dn/d\nk = for pre-WMAP datasets (see text, dotted 
contours), WMAP three-year data (dashed contours) WMAP five-year data (soHd contours), and WMAP five-year data plus 
ACBAR (filled contours). The pre-WMAP limits were generated using the same parameters as the WMAP contours, except 
that the reionization optical depth was constrained to be r < 0.2. 



